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ABSTRACT 

The compact Abelian model in three space-time dimensions is 
studied in the presence of external electromagnetic fields at finite 
temperatures. We show that the deconfinement phase transition is 
independent on the strength of the external fields. This result is in 
agreement with our observation that the external fields create small- 
size magnetic dipoles from the vacuum which do not infiuence the 
confining properties of the model. Contrary to the deconfinement 
phase, the internal field in the direction of the applied external field is 
attenuated in the confinement phase, this screening becomes stronger 
with decreasing temperature. 
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1 Introduction 



Compact Abelian gauge theory in three Euchdean dimensions is proven to possess the 
property of permanent confinement 0, due to presence of Abehan monopoles: any pair 
of test electric charge and anti-charge are confined by a hnear potential. The monopoles 
are topological defects which appear due to the compactness of the gauge group. In three 
dimensions the monopoles are instanton-like. 

The confining property of the model is lost at sufficiently high temperature. The 
confinement-deconfinement phase transition — which is expected ^, |^ to be of Koster- 
litz-Thouless type p — was studied on the lattice both analytically ^ and numeri- 
cally 1^. A thorough numerical analysis of the phase transition was done in Ref. 0] 
where it has been demonstrated that the monopoles are sensitive to the transition. In 
the confinement phase the monopoles are observed in the plasma state while in the de- 
confinement phase the monopoles appear in the form of a dilute gas of magnetic dipoles. 
Similarly to the monopole plasma the dipole vacuum, although not confining, still has a 
non-perturbative nature In the confinement phase both monopole density and string 
tension differ from semiclassical estimates []I| which neglect monopole binding. However, 
an analysis of monopole clusters shows that the relation between the string tension and 
the density of monopoles in magnetically charged clusters is in reasonable agreement with 
those predictions. 

An alternative explanation of the deconfinement phase transition, based on Svetit- 
sky-Yaffe universality arguments |^, was given in Refs. ^. The phase transition was 
demonstrated to be accompanied by restructuring of the U{1) vortex system in an effective 
2D Abelian spin model. In the confinement (low temperature) phase the vortices exist in 
the plasma state while in the deconfinement (high temperature) phase the vortices and 
anti-vortices form bound states. Both monopole and vortex |jT^ binding mechanisms 
have been studied in the context of the Georgi-Glashow model, too, a limiting case of 
which is the compact Abelian gauge model. 

In this paper we continue to consider the deconfining mechanism by monopole pairing 
which seems to have interesting counterparts in realistic gauge theories. The formation of 
monopole pairs is qualitatively similar to the binding of instantons in instanton molecules 
with increasing temperature in QCD suggested to be responsible for chiral symmetry 



restoration |TT|. An external magnetic field affects the phase diagram of the non- Abelian 
theory as was shown both analytically |T^ and numerically |13[ . In the electroweak theory, 
the formation of Nambu monopole — anti-monopole pairs, a remnant from a dense medium 
of disordered Z-vortices and Nambu monopoles which characterizes the high temperature 
phase, is accompanying the transition towards the low temperature phase [Q. The effects 
of the external fields on the phase transition temperature and the electroweak baryogenesis 
in the electroweak theory were discussed in Ref. [|l^ . 

In three dimensional gauge theories the inclusion of external fields allows to study the 
vacuum energy density in the SU{2) gauge theory [1^, vortex dynamics in the Abelian 
Higgs model [O] and dynamically generated fermion mass in the Abelian gauge theory 



with fermionic |]18| fields. Here we study the infiuence of the external electromagnetic 
field on the confining and monopole properties of the compact Abelian gauge model in 
three dimensions. 

The plan of the paper is as follows. In Section || we describe the basics of the lattice 
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formulation of the compact U{1) model in 2 + 1 dimensions and how the external field is 
implemented. The screening of the external electric and magnetic fields and the properties 
of the Polyakov loop correlators in the presence of those fields are studies in Section |^. 
Differences and similarities between magnetic and electric fields are stressed. The string 
tension in the presence of the external fields with a flux next is investigated in Section ^ 
and the phase diagrams in the next — P plane are presented in Section |, separately 
for magnetic and electric fields. In the case of an applied electric field we see that the 
so-called "bulk Polyakov loop" cannot play any longer the usual role in localizing the 
deconfining transition. Results for various monopole properties are collected in Section |^. 
Our conclusions are summarized in the last Section. 



2 The compact lattice model in external fields 

We study 2 + 1 dimensional compact lattice electrodynamics in constant external electro- 
magnetic fields. Following Ref. []T9[ we use the action 



where 6p is the U{1) field strength tensor represented by the curl of the compact link 
field 9i, and is the field strength corresponding to the external field. P is the lattice 
coupling constant related to the lattice spacing a and the continuum coupling constant 
of the 3D theory as follows: 

P = —,- (2) 
a 93 

At finite temperature the lattice is asymmetric, x Lt, Lt < Lg] Li = L2 = Lg and 
L3 = Lt are the spatial and temporal extensions of the lattice, respectively. In the limit 
Ls —>■ 00 the temporal extension of the lattice Lt is related to the physical temperature, 
Lt = l/{Ta). Using eq. (|^) the temperature is given via the lattice parameters as follows: 



L- H 
gi Lt 



2-r ■ (3) 



Note that in 2 + 1 dimensions there is no symmetry anymore between the three com- 
ponents of the field strength tensor. The closest relative of the true magnetic field is 
F12 distinct from the others, while there is still a symmetry between F13 and F23. With 
this distinction in mind one can conditionally call them the "magnetic" and "electric" 
components of the field strength tensor, respectively. As long as one does not introduce 
external fields, even at finite temperature there was no need to distinguish between them. 

Without loss of generality we represent an external electric field as a field directly 
coupled via (|l]) to plaquettes lying in the 31 plane. The value of the field is quantized ||T9| : 



E = eif = p^, ueez, 9'if = e^f = o. (4) 
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To study the influence of a magnetic field (and to draw parallels to analogous studies in 
AD [^) we implement a magnetic field 5 as a field directly coupled via (jl|) to plaquettes 



lying in the 12 plane, 



B = eif = '^, umez, e^f = e^f = o. (5) 

The quantization of the external fields (^), (|^) is necessary to match with the periodic 
boundary conditions imposed on the lattice |T3 . Indeed, these conditions imply that the 



rectangular Wilson loop of the size Li x Lj in the ij plane must be equal to unity. On the 
other hand the Wilson loop is equal to exp{i$}, where $ is the flux going through the ij 
plane. Therefore, we get the quantization condition for the flux, $ = 27rnext, '^ext £ 
which implies in turn the relations (|^), (^) for the constant electric or magnetic field, 
respectively. We use the notation he/m instead of next where we want to emphasize the 
electric/magnetic nature of the external field. 

Naturally, from the form of the action (0) the number of external fluxes is restricted, 
< n^xt < LiLj/2. In our analysis we restrict that number to < ricxt < LiLj/A. One 
should note however, that the largest considered fluxes correspond to lattice artifacts. 
Indeed, the physical electromagnetic field must have a strength smaller than the scale 
corresponding to the lattice ultraviolet cutoff, a~^. This leads to the restriction of the 
lattice field strength: 6'cxt ^ 1? or, next ^ ^max = [LiLj/{27i)], to avoid too large external 
fields. 

As in the case of zero field we restrict ourselves to a finite temperature lattice 32^ x 8 
varying the strength of the constant external electric or magnetic field and considering 
the lattice gauge coupling range 0.1 < /5 < 3. For this particular lattice size we get the 
upper bounds for the electric and magnetic fields: n^^"^ = 40 and n™/"" = 163, respectively. 
From our studies at zero external fields we |^ found the deconfinement phase transition 
at (3c = 2.346(2) using the Polyakov loop susceptibility in satisfactory agreement with 
Ref. i. 

We briefly recall the Monte Carlo algorithm used and described already in Ref. ^ . The 
algorithm combines a local Monte Carlo step with a global refreshment step to improve 
ergodicity. The local Monte Carlo algorithm is based on a 5-hit Metropolis update sweep 
followed by a microcanonical sweep. Following Ref. |Jl9[, the global refreshment step 
consists of an attempt to add ari additional unit of flux of randomly chosen direction and 
sign {i.e. adding a gauge field 9i constructed over the whole lattice to the previous one, 
l^i + ^i]mod27r' subjcct to a global Metropolis acceptance check. The acceptance rate 
for a global update does not depend on the external flux next and it is roughly equal for 
global shifts of the electric or magnetic flux. However, it changes with (3, reducing from 
0.58 ... 0.6 at /3 = 1.0 to 0.18 . . . 0.20 at /3 = 2.9. 

3 Polyakov loop correlators and screening 

Usually, the Polyakov loop 

Lt 

L(x) = expji J dt Ao(x,t)} (6) 
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is used as a the basic quantity to probe the confining properties. The integration runs 
along a loop parallel to the time axis and located at the spatial 2D coordinate x. The 
Polyakov loop operator inserts an infinitely heavy test particle with unit electric charge 
into the vacuum of the theory. The v.e.v. of this operator expresses the free energy 
F of the inserted particle, (L(R)) = exp(— F/T) and is an order parameter to signal 
deconfinement. 

In the absence of an external field the correlation function between two Polyakov 
loops can be expressed via the interaction potential V{R) between infinitely heavy electric 
charge and anti-charge: 

(L(0)L*(R)) (X exp{-Lt V{R)} . (7) 

The leading behavior of the potential V^(R) in the low temperature phase corresponds to 
a linearly rising potential, V{R) = a R, where a is the tension of the string between the 
test electric charges separated by the distance R = |R|. In the high temperature phase 
the potential is of Coulomb type (it rises logarithmically with increasing distance between 
test particles). 

Now we discuss the Polyakov loop correlations in both electric and magnetic external 
fields. We introduce a 3D vector combining electric and magnetic components by 

F = {F23, F31, F12) = {Ex, Ey, Bz) . (8) 

As we will immediately see, an external (constant) electric field substantially modifies the 
behavior of the Polyakov loop correlator (^ already on the tree level. 

Indeed, the electromagnetic field inside the medium can be decomposed into two parts, 
F™* + F'^. Here F™* is the mean field inside the medium, which is non-zero due to the 
presence of the external electromagnetic field. The field F^ corresponds to the quantum 
fiuctuations around the mean value F™*. Due to the Abelian nature of the Polyakov loop 
the correlation function (|^ can be written as product of two contributions 

(L(0)L*(R))^.., = (exp{z/dx^A^(x)})^ = e**'^^^") ■ G(R, F^^*) , (9) 
G(R,F-*) = (L(0)L*(R))l,, = (e^*^(^^'));^..., (10) 
where C = Iq — 1-r, is the contour corresponding to the test particle trajectories, and 

^c{F) = / dx,A,{x) = I da^F^{x) , (11) 

is the fiux of the electromagnetic field F which goes through the surface Sc spanned on 
the contour C. The subscript F"^^^ in eq. ( p!oD indicates, that the vacuum expectation 
value of the quantum part of the Polyakov loop has been taken with the action eq. (|l|) 
corresponding to non-zero external fiux. 

The correlator (p!OD is taken over quantum fiuctuations only, as indicated by the super- 
script q. Note that by definition the field A'^ fiuctuates near the minimum of the action. 
Thus negative modes are absent and the potential V^(R; F^"^^) oc — logG(R, F^^^)/Lt must 
be real. The dynamics of the monopoles affects the quantum fiuctuations and consequently 
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the potential V. The external fields disturb the monopole medium and therefore we may 
expect that in general the inter-particle potential V may get infiuenced by the external 
fiuxes. 

Summarizing, the quantum average of the Polyakov loop in the presence of the external 
constant electromagnetic field is given by the following formula: 



(L(0)L*(R))^.., oc expj^^clF''^*) - L,\/(R;F-*)} , 



(12) 



where a non-vanishing Rnx pierces the surface spanned by the contour C and 

the trajectories of the test particles are placed along the (temporal) ^-direction. 

The properties of the internal fields inside the system are important to understand 
the behavior of the Polyakov loop correlators. For external fields given by eqs. (^), (^) 
the total external fluxes through the corresponding planes are quantized. Analogously, 
the internal electric or magnetic fluxes through any plane in any gauge fleld conflguration 
are quantized as well (the considerations similar to those in Section ^ are applicable in 
this case as well). Note that the global update step is consistent with the internal flux 
quantization since the update allows to change the total flux of the system by just one 
unit. 

From now on we study only two possibilities for the external fleld, an electric fleld in 
y direction and magnetic fleld in z direction. 



(0, E, 0) , or F^"* = (0, 0, B) 



(13) 



It is clear that only the electric component contributes to the flux $c(F™*) in eq. (|T2|). 
In our case the mean internal electric and magnetic flelds which are the actual average 
flelds present in the medium are given by {O^i) and (6*12), respectively, where 



(14) 



s f X 



Following Ref . |21| we show in Figure |I] the normalized average internal electromagnetic 
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Figure 1: The normalized internal (a) electric and (b) magnetic flelds vs. j3 at flxed next- 
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fields (^3i)/^ = F'^^/E and {612) / B = F'^^ / B along the directions of the applied fields as 
function of the gauge coupling (3 for various external field fiuxes next- One can see that in 
the confinement phase the internal fields are much weaker than the external fields while in 
the deconfinement phase the internal and external fields almost coincide with each other. 
The measured average internal fields in the directions transverse to the applied external 
field remain zero. 

The attenuation of the external fields in the confinement phase resulting in weaker 
average fields in the medium happens due to the monopole plasma [||] and can be called 
screening: the monopoles produce a finite correlation length ^. If we would apply an 
external field to the box with the monopole plasma using frozen boundary conditions 
then the field would be greatly suppressed inside the media at distances (from the box 
boundaries) larger than S,- Thus the averaged field inside the box with the plasma should 
be smaller than the external field. The smaller the correlation length the smaller the 
averaged field inside the media should be. 

As the temperature increases the monopoles form more and more magnetic dipole 
bound states. In the dipole plasma the screening is absent and the correlation length 



is infinite. The internal field is expected to be equal to the external one what is clearly seen 
in Figures |T| for external fields of small strength. Moreover, the screening of the external 
magnetic and electric field with the same strength (related to each other as n^; = hm/ 
4 due to difference in the lattice extensions in spatial and temporal directions) is very 
similar for weak external fields. This happens because the nature of the monopoles is 
not purely "magnetic" at finite temperature. At large magnetic fields the response of the 
media is different in magnetic and electric cases. In Section ^ we show that the observed 
features of the screening are qualitatively related to specific properties of the monopoles 
and their bound states. 

Let us consider the properties of the correlator ([T2|) on a lattice of finite size for a non- 
vanishing external electric field E using the lattice Polyakov loop at position x = (x, y) 

L(x) =exp(^zX]^3(x,^)^ . (15) 

We put the test charge and anti-charge at the points (0,0) and {x,y), respectively. The 
fiuxes $c(^3i) going through the surface spanned between the test particle trajectories 
depend on the internal electric field ^31, they are quantized and peaked around an average 
value. Therefore, they can adequately be described by taking into account only the 
"most probable" fiux state labeled by the integer njnt = nint(n,ext, -^i) which depends 
on the strength of the external field, temperature and lattice geometry. In that case 
$c = 27rninta; /Lg and -Emt = (^^31) = STrriint /{LgLt). The Polyakov loop correlator reads 
as follows, 

(L(0, 0) L*{x, y))j, oc exp{27rzni,t ^ ^ Lt V{x, y-E)], (16) 

with the oscillating part of the correlator defined by the electric component of the internal 
field. 

Eq. (|16D simultaneously characterizes both the screening of the external field (phase 
factor) and the potential of the test electric charges (modulus). We expect that in the 
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deconfinement phase the external field is not screened and the correlator is dominated by 
^int ^ ue- In the confinement phase the field must be screened and nint becomes smaller 
with decreasing /? (temperature). 

Thus the Polyakov correlator in the deconfinement phase is an oscillating sign-chan- 
ging quantity at fixed y coordinate as function of the x coordinate. The maxima of the 
oscillations become damped with increasing x due to the Coulombic interaction between 
the test particles. 

The correlations at fixed x and changing y should be of the same sign and decay with 
increasing y. This behavior is nicely illustrated by the (real part of) the Polyakov loop 
correlators (|16D in the x — y plane shown in Figure 0(a) for the deconfinement phase. The 




(a) (b) 

Figure 2: Real part of the Polyakov loop correlator in the x — y plane ([T6|) for an external 
electric fiux n^; = 4 in (a) deconfinement, (3 = 2.6 and (b) confinement, (3 = 2.0, phases. 
Half of the lattice 32^ x 8 is shown. 

measurement was performed in an external electric field with fiux n^; = 4. 

In the confinement phase we expect that the sign fiuctuations of the correlator should 
be strongly suppressed due to area-law decay. Thus the correlator must go to zero rapidly 
with increasing distance between charge and anti-charge, a/x^ + y^. The corresponding 
correlator, presented in Figure ^(b) at /5 = 2.0, shows the expected suppression. In 
addition, we observe a small local maximum at x ~ 8 and y = arising from a fiux of 
value nint close to 4 in agreement with the behavior of the internal fields in Figure Q for 
this relatively large f3 value. 

Finally, let us comment here that an induced internal magnetic field (parallel to its 
parental external field) which is directed along the Polyakov loop cannot contribute to 
the "classical" phase $c cind, therefore, oscillations cannot appear. Thus, in an external 
magnetic field, the correlator is rapidly approaching zero with increasing spatial distance. 
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4 String tension in the presence of fields 



According to Ref. the point-point correlation of the two Polyakov loops in the absence 
of an external field can be parametrized as follows^: 

Qipi(xi-x2)+ip2{yi-y2) 

{L{xi, yi)L*{x2, 2/2)) = const ■ ^ ■ , . + . . . , (17) 

^ 1 - cospi - C0SP2 + cosh{aLt) 

where the sum runs over all possible momenta, pi^2 = 0, . . . , 27r(Ls — 1)/Ls and a denotes 
the "temporal" string tension. 

To evaluate the string tension, we use two Polyakov "plane" operators which are 
defined as a sum over Polyakov loop positions along and perpendicular to the external 
field. As an example we take the electric field in y direction. Then 

L|,(x) = X;L(x,i/), L^{y) = Y^L{x,y). (18) 

y=l x=l 

The correlator of the plane-plane correlators may be written as a sum over the corre- 
sponding positions of the point-point correlation functions, 

(L||(0)L|(x)) = {my^)L*ix,y2)), (19) 

1/1,2=1 

(L^(O)Ll(y)) = {L{x^,0)L*{x2,y)). (20) 

a;i,2=l 

In the presence of an electric field E along the y axis this correlator must be modified 
in accordance with our discussion in Section || 



{L{xi,yi)L*{x2,y2)) {L{xi,yi)L*{x2,y2)) 



E 



Xi - X2 



exp ( -27rmint — ) {L{xi, yi)L*{x2, y2))„^ • (21) 



The subscript ue indicates that the Polyakov loops are calculated in the state with an 
external electric flux n^; in the y direction. The internal electric flux riint depends implicitly 
on that external flux. Combining eqs. ([T7|), (|2T| ) with eqs. (pJf), ( pOD and taking the sums 
over the momenta pi^2 explicitly, we get: 



:22i 



(L||(0)L|(a;))^^ = const ■ e^™-'^/^- coshJfxLt (a; - y 

{L^{0)L*_^{y))^^ = const ■ cosh acsicr,nint) Lt(y - y) ' (23) 

where formally a string tension coefficient crefr(cr, nint) for the plane-plane correlator per- 
pendicular to the fleld is: 

cTefr(o", nint) = -— arccosh cosh(crLt) — cos(27raint / Lg) + 1 . (24) 



L 



''Here and below we take into account only the lowest "mass state" corresponding to the actual string 
tension tension, unless specified otherwise. 
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Thus in the presence of an external electric field the plane-plane Polyakov loop correlator 
parallel to the electric field oscillates with a decreasing amplitude. 

The plane-plane correlator perpendicular to the field decreases exponentially (without 
oscillations) as function of the distance between the planes. The decrease of this correlator 
is controlled by an effective string tension (|^) which is a function of the external electric 
field (via nint) and the actual string tension. The essential message here is that this 
effective string tension does not tell anything about confinement properties, since the 
confinement is described by the real string tension a. Note that the effect of the external 
electric field is absent if the internal fiux is equal to quantized values 



n 



int 



N e 



(25) 



The origin of the a^s dependence on the external electric field strength is related to 
the oscillatory behavior of the Polyakov loop correlator (|T6|). Indeed, the bulk correlator 
{\L\'^)e (defined below in eq. (p6D) comprises all possible Polyakov loop correlators. At 
small and increasing nint the sum of the oscillating quantities leads to a suppression of 
{\L\'^)e while at riint satisfying the special condition ( ^3]) the oscillations disappear and as 
a result the effective string tension ( |2^ coincides with the string tension in the absence 
of the external electric field. 

To get rid of the direct infiuence of the electric field on the string tension ( P^ we 
use only the real part of the Polyakov plane-plane correlator which is parallel to the 
external field. We fit the numerical data by the real part of eq. ( ^2]) using (besides trivial 
prefactors) a and riint G ^ as fitting parameters. At the considered values of /3 > 1.6 we 
find that the best fit is given by nint = ""-ext in accordance with our discussion in Section ||. 

In the absence of external fields, the dependence of the string tension on /3 has been 
found in Ref. and it is shown in Figure ^(a). The string tension is a rapidly decreasing 
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Figure 3: (a) The string tension vs. (3 without external field taken from Ref. 0. (b) 
Fitted string tension for various (3 values as function of the external electric fiux ue- 



function: it is non-zero in the low temperature phase (low (3) and becomes very small (due 
to finite size effects) in the high temperature phase. From Figure |](b) we conclude that 
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within errors the string tension does not depend on the external flux ue- This observation 
allows us to conclude that an external electric field in the compact {2 + l)D Abelian gauge 
theory does not change the confinement behavior compared to the zero field case. 

For an external magnetic field (in the "temporal" z direction) the oscillating factor in 
the Polyakov loop correlator is absent. Therefore, we have used eq. (p2D with n^^t = 
to fit the plane-plane correlation functions of the Polyakov lines. The numerical results 
for the string tension a are similar to those in the case of the electric field shown in 
Figure 0(b): within errors of the fit values we do not observe a dependence of the string 
tension on the strength of the external field. Therefore we do not show the plot of the 
string tension vs. the external magnetic fiux. 

Summarizing, both magnetic and electric fields do not infiuence the string tension. 
This fact is analyzed in terms of the monopole properties in Section ^. Before doing so 
we discuss the phase structure of the model using the Polyakov loop expectation value 
and the corresponding susceptibility. 



5 Phase structure 

The confinement-deconfinement phase transition is usually detected using the Polyakov 
loop vacuum expectation value. It is convenient to study f.e.f.'s of the powers of the 
"bulk Polyakov loop" defined as follows: 

(1^1) = (I E^WI) ' (1^1') = :^ (I E^WI') = (E^W^^(y)) • (26) 

s X ■'^s X ■'^s x,y 

The Polyakov loop susceptibility is expressed via these quantities: 

XL = {\L\f - {\L\') . (27) 

According to the free energy arguments in the deconfinement phase the quantity \L\ 
should be non-vanishing while in the confinement phase this quantity becomes small (it 
approaches zero in the infinite volume limit). 

In the following we show that the expectation value of the Polyakov loop, similar to 
the correlator discussed above, also gets a large classical contribution due to the external 
electric field. Therefore it should not be used in general as an order parameter to probe 
the restoration of confinement (except for some special values of the electric field). 

5.1 Electric field 

The behavior of the Polyakov loop vs. (3 aX various values of the external electric field 
is shown in Figure ^(a). The plot of the Polyakov loop at zero field agrees with general 
expectations. However, as the electric field is turned on, the vacuum expectation value 
of the Polyakov loop is decreasing. One might conclude that the external electric field 
restores the confinement phase (while being in deconfinement at n^; = 0) which is in clear 
contradiction with the results for the string tension shown in the previous Section. We 
remind the reader that the external electric field modifies the properties of the Polyakov 
loop correlator classically. In particular, for not too large it enhances what we called the 
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Figure 4: (a) The absolute value of the bulk Polyakov loop p6| ) and (b) the Polyakov 
loop susceptibility (|27|) vs. [3 for various values of the external electric field, he- 



effective string tension measured from the Polyakov plane-plane correlators perpendicular 
to the electric field (c/. eq. (|2^)). But, as we noted above, this effective string tension 



does not describe the confining properties of the system. 

Thus we may conclude that the rapid decrease of the Polyakov loop as function of 
He also does not mean that confinement is restored. Note that at large special values 
of the internal flux nint, eq. the module of the Polyakov loop and the susceptibility 
(Figure ^(b)) as function of (3 do not differ from the zero field case. 

Let us estimate the classical correction to the Polyakov loop expectation value due to 
the external electric field E. For this purpose we consider the squared modulus of the 



Polyakov line, eq. (|26|). In the presence of E it can be written as follows: 



(|LH^ = const . ^e^^^^"-^/^^ (L(0, 0)L*(x, y)) . (2^ 



Taking into account eq. ([TTp for the point-point Polyakov loop correlator and summing 
over all momenta pi in this equation we getQ: 

{\l?)e=T.^, — 77% — \;r-r^- ^^9) 

^0 cos(27rnint/lvs) - cosh(am-t^t) 

Here the expansion over the excited mass states cr^ is written explicitly and the lowest 
state corresponding to the string tension: (Tq = a. From Figures § it follows that xl 

for all 13 and ue values. Thus the square root of eq. (p9D can serve as a good 
estimator for the Polyakov loop quantum average {\L\). 

We use only the first two terms of the expansion (pOf) to fit the behavior of the measured 
Polyakov loop vs. the external electric flux number. The corresponding fits for the 
confinement and deconfinement phases are shown in Figure ^(a,b), respectively. The best 

^ According to our considerations above the dominating internal flux is equal to the external flux in 
the vicinity of the phase transition and thus we can safely put riint = ns- 
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Figure 5: The absolute value of the bulk Polyakov loop, p6| ) vs. the external electric 
flux He and its fit by formula (^91) in (a) the confinement phase (/3 = 1.6) and (b) the 
deconfinement {(3 = 2.5) phase. 



fits for the string tension in the confinement phase (at f3 = 1.6) give: do = 0.167(13) 
and (Ti = 0.85(6). The value of ctq can be compared with the string tension without 
external field, a = 0.170(4). These values agree with each other within statistical errors. 
In the deconfinement phase (/3 = 2.5) the same picture is found: ctq = 0.0096(5) and 
0"! = 0.18(1), while the independent measurements at zero field lead to cr = 0.0101(1). 
Therefore, the behavior of the Polyakov loop approximated via eq. (|29| ) in the external 



electric field is consistent with the string tension measurements at zero-field. 

So, in the case of non-vanishing external electric field, the bulk Polyakov loop expec- 
tation value may vanish regardless of the value of the actual string tension. Indeed, even 
at cr = (which is not accessible on the finite lattice due to finite volume effects) the 
squared Polyakov loop expectation value (pOf ) is decreasing when turning on the external 
electric flux. As we have mentioned in Section ^ the influence of the external field on 
the Polyakov loop correlations is absent provided the condition (^) for the internal flux 
nint is fulfilled. Now, a similar effect is observed for the Polyakov loops themselves: the 
corresponding values for ue = nint = 32, 64 clearly coincide with he = ^ data. Thus 
we conclude that the observables based on the Polyakov loops may have a usual physical 
sense only in the cases of the quantized internal field (pS]). 

Analogously, Figures §(a,b) show that the dependence of the Polyakov loop and its 
susceptibility on the coupling constant (3 are the same (within errors) for n^; = , 32, 64. 
The peaks in susceptibility may serve as good indicator of the (pseudo) critical coupling 
constant. At these distinguished values of the external electric field we have fitted the 
susceptibility near its maximum by the following function: 

^■^'f^-^fTFTTy.. (30) 

where q and j3c are fitting parameters. The results are shown in Figure ^(a) in the ue — (3 
plane. One can see that the influence of the external field on the critical temperature is 
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Figure 6: Phase diagrams in the ue/m — P plane derived from the module of the Polyakov 
loop for (a) electric and (b) magnetic fields, respectively. 



negligible. 



5.2 Magnetic field 

The tree level contribution to the Polyakov loop observables is absent in the case of the 
external magnetic field. This is confirmed by Figures 0(a,b) where the Polyakov loop 
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Figure 7: (a) The Polyakov loop and (b) its susceptibility at various values of the external 
magnetic field um vs. p. 



and its susceptibility are shown as functions of (3, respectively. The data points for all 
considered values of the external magnetic flux coincide with each other within errors. 
We fit the Polyakov loop susceptibility by eq. (0) to get the (pseudo) critical couplings. 
The corresponding phase diagram is shown in Figure |](b) in the um — P plane. One sees 
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that the phase transition points are insensitive to the strength of the external magnetic 
field. 

It is interesting to note that in non-Abelian gauge theories in 3 + 1 dimensions the 



situation is quite different [0. Here the position of the deconfining phase transition 
depends strongly on the strength of the applied magnetic field. The simplest, somewhat 
naive explanation of this fact could be as follows. In the non-Abelian gauge theory the 
coupling between gauge fields is stronger than in the Abelian theory. In particular, there 
is a correlation between spatial and temporal components of the fields which may imply 
that the external magnetic field (encoded in the spatial components) induces internal 
electric fiuxes on the quantum level. As we have seen the electric fiuxes strongly infiuence 



the Polyakov loop correlators generating the effective string tension (|24D. Contrary to 
the Abelian case this influence is not a merely classical (or, just inherent to the way of 
introducing the external fleld) but a real quantum effect. 

In the next Section we study the effects of the external flelds on the agents of conflne- 
ment, the Abelian monopoles, to explain the observed behavior of the system. A similar 
study for the non-Abelian gauge theory is underway . 



6 Monopole properties 

The basic quantity to describe the behavior of the monopoles is the monopole density, 

Pmon = , (31) 

c 

where rric is the integer valued monopole charge inside the cube c deflned in the standard 
way P5[ : 



Pedc 

In our previous study |0 we have demonstrated that the monopoles are sensitive to the 
phase transition in the compact Abelian gauge model at flnite temperature. Although 
in (2 + 1)D we have magnetic and electric fleld among the three components of the fleld 
strength tensor, the sources of the respective fluxes will be simply called "monopoles" or 
"magnetic charges" in the following. 

The mechanism which drives the flnite temperature deconflnement phase transition 
is monopole binding. In the zero temperature case the plasma of monopoles and anti- 
monopoles can explain the permanent conflnement of oppositely charged electric test 
charges |jl| in bound states, kept together by a linear potential. Conflnement appears due 
to the screening of the magnetic fleld induced by the electric current circulating along the 
Wilson loop. Monopoles and anti-monopoles form a polarized sheet of flnite thickness 
("string") along the minimal surface 1^41 spanned by the Wilson loop. The formation of 
the string leads, for non-vanishing electric current, to an excess of the free energy equal 
to a\A\. 

At flnite temperature, dipoles are formed both in the conflnement and deconflnement 
phases. In the deconflnement phase tightly bound dipoles dominate in the vacuum. The 
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dipole plasma is inefficient to completely screen the field created by the electric currents 
running along the pair of Polyakov loops. This explains the absence of confinement in 
this phase. 

Besides measuring the density of monopoles, we have studied the properties of the 
monopole ensembles by investigating the structure of monopole clusters. Clusters are 
connected groups of monopoles and anti-monopoles, where each object is separated from 
at least one neighbor belonging to the same cluster by a distance less or equal than -Rmax- 
In the following we use -Rmax = ct which means that neighboring monopole cubes 
should share at least one single corner[|. The increase of the coupling constant leads not 
only to an increase of the temperature, eq. (0), but also to a decreasing lattice spacing a, 
eq. (0). Thus at different (3 the same characteristic distance i?max corresponds to different 
physical scales. Therefore our results presented here are of a qualitative nature. 

In our study without external field we have found that the dipoles are oriented dom- 
inantly in the temporal direction. At the confinement phase transition mostly clusters 
with two constituents or single monopoles and anti-monopoles were observed. Decreasing 
further the temperature (or (3), the monopoles become dense and start to form connected 
clusters (on a coarser and coarser lattice) containing various numbers of monopoles and 
anti-monopoles. The largest clusters have been found to be more and more spherical. 

Finally we observed that only charged monopoles clusters in the plasma (mainly indi- 
vidual monopoles) are needed to explain the measured string tension and, therefore, are 
responsible for confinement. 

6.1 Electric field 

To begin, we plot in Figure |](a) the total monopole density p as function of the coupling 
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Figure 8: Monopole density pmon (a) vs. [3 for various external electric fluxes, n^; (b) vs. 
ue for various [3 values. 



P for various values of the external electric flux he- The monopole density is a decreasing 

■^In Ref. |p6|| a similar definition has been used to investigate tightly packed clusters with Rmax = a- 
In our case the condition for the cluster is more relaxed. 
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function of (3 at any value of the external field. However, p increases as the function of 
the strength of the applied external field, Figure ||(b). The effect of the external field is 
very essential: the monopole density is increased up to almost two orders of magnitude 
(depending on the temperature) for the largest external flux values (compared to the 
system at zero external field). 

However, an increased total monopole density does not mean in general that the confin- 
ing properties of the system are enhanced. Only charged monopole clusters in the plasma 
state contribute to the string tension between electrically charged test particles. Tightly 
bound monopole pairs are not expected to contribute to the string tension. Therefore, we 
use our cluster labeling algorithm to look into the structure of the monopole ensemble. 
In Figures Ka,b) we show how the densities of single monopoles (clusters made of just 
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Figure 9: The density of clusters of size (a) = 1 (single monopoles and anti-monopoles) 
and (b) iV = 2 (dipoles) vs. external electric flux ue at various values of (3. 



one (anti-)monopole) and of dipoles (clusters made of two oppositely charged objects), 
respectively, depend on the strength of the external electric field. 

One clearly sees that the plasma component of the single monopole ensemble does 
not feel the electric field at all. Thus the confining properties of the system should 
not depend on the external electric field in agreement with the conclusions made in the 
previous Sections. 

On the other hand, the dipole density changes drastically with increasing external 
field: the field creates the magnetic dipoles from the vacuum. Note that the larger the 
temperature (or, equivalently, /?), the larger is the increase of the dipole density. This 
fact is connected with the screening of the external fields inside the medium discussed in 
Section |^. The larger the temperature the larger the fields are inside the medium. As a 
result, the effect of the external field becomes stronger with increasing temperature. 

In a non-zero electric field the system is anisotropic in all directions. The electric field 
is directed along the y axis, the "temperature" direction, 2, is influenced by compacti- 
fication. Therefore, we have to expect that the average sizes of the dipoles in different 
directions are not the same. At zero or small external field the dipoles are mainly directed 
along the temporal axis 0. 
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Increasing the external field, the dipoles are expected to become elongated along the 
direction of the applied field. Moreover, we have observed in Section ^ that the strength 
of the internal field inside the medium relative to the external field increases as function 
of the coupling (3. Thus the elongation of the dipoles in the field direction should increase 
with (3. 

All these effects are demonstrated in Figure using eUipsoids, the axes of which are 




(a) (b) 

Figure 10: The mean dipole anisotropy for increasing values of (a) the external electric 
flux He at fixed P = 1.8; (b) (3 at fixed ue- 

equal to the average dipole sizes in the x, y and z directions. In Figure p!0|(a) we show how 
the mean dipole anisotropy changes with increasing external electric flux at fixed (3 = 1.8 
(in the confinement phase). Figure |l3(b) demonstrates the dependence of the ellipsoids 
on j3 at fixed flux ue- For convenience the projections of the ellipsoids onto the x — y and 
X — z planes are also shown. 

6.2 Magnetic field 

The influence of the magnetic field on the monopole densities is very similar to that of 
the electric field. In Figures |ll](a,b) we present the measured total monopole density and 
the dipole density, as functions of f3 and the external fluxes, respectively. To compare 
the measurements for magnetic fluxes with those at non-zero external electric field, we 
note that for our lattice geometry the strengths of the magnetic and electric fields are 
equal to each other provided hm = '^'tie- According to those Figures both total monopole 
density and dipole density do not depend on the type of the external fields if the field 
strengths are the same. We have checked that the single monopole density coincides for 
magnetic and electric fields of the same strength as well. Thus the cluster structure does 
not depend on whether the external field is of electric or magnetic kind. 
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Figure 11: (a) Total monopole density pmon vs. (3 for various external flux values next- 
(b) density of dipole clusters vs. ricxt = = um/^: for some /3 values. 



The mean dipole anisotropy depends on the type (direction) of the external fleld. We 
show the behavior of this quantity in Figures 0(a,b) for increasing external magnetic 




(a) (b) 

Figure 12: The mean dipole anisotropy for increasing values of (a) the external magnetic 
flux hm at flxed (3 = 1.8; (b) j3 at flxed um- 

fleld and (3, respectively. In the external magnetic fleld the system is symmetric in the x 
and y directions. The dipoles become more elongated in the (temporal) z direction for 
stronger external magnetic flelds. Increasing the coupling /3 leads to a larger anisotropy 
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since the medium does not screen the external field in the deconfinement (large/3). 

The orientation properties of the dipoles infiuence the resulting electromagnetic fields 
in the medium, Figure |l|. According to Figures |l3(a) and [l^(a) the stronger the external 
field the larger elongation of the dipoles in the direction of the field is. In turn this 
increases the dipole momenta in the direction of the applied field and, as a result, leads to 
an enhanced screening of the external fields by the media. This effect is clearly observed 
in Figure 0. 

Another interesting effect due to the dipole orientation is the clear difference in screen- 
ing for strong magnetic and electric fields of equal strengths in the deconfinement phase. 
Figure |I]. According to Figure 0(b), the dipole densities for both cases are the same. 
Since Lt < Ls, the dipole magnetic moment density projected to the magnetic field direc- 
tion is larger than the corresponding quantity for the electric field direction. On the other 
hand, the larger the dipole moment density inside the medium the stronger the field 
attenuation is. Thus at strong external fields the screening must be more effective for 
magnetic fields compared to electric fields. This mechanism works at high temperatures 
(large (3) where the dipole fraction is dominant (c/. Figures ^(a,b)). 

7 Conclusion 

We have investigated the properties of the 3D compact electrodynamics at finite tem- 
perature in external constant electric and magnetic fields. The main result is that the 
deconfinement temperature is insensitive to both electric and magnetic external fields. We 
have found the reason for this behavior in terms of the monopole degrees of freedom: the 
external fields create tightly bound magnetic dipoles (monopole — anti-monopole bound 
states) from the vacuum while the density of unpaired monopoles (which are responsi- 
ble for the confinement of electric charges) stays unchanged. This result is not obvious 
from the beginning since another option is possible: the external fields could destroy the 
monopole bound states enforcing confining properties of the medium. This is not the 
case. 

At zero external fields the magnetic dipole states are more elongated in the tempera- 
ture direction. The external magnetic field which is parallel to the temperature direction 
makes this elongation stronger. However, an external electric field turns the dipoles to 
the corresponding spatial direction. The effects of the external field on the medium are 
stronger in the deconfinement phase in which both electric and magnetic external fields 
are not screened. 

We have also shown that the external electric field infiuences the Polyakov loop clas- 
sically (or, in other words, on tree level). This leads to a vanishing Polyakov loop and, in 
certain cases, to a non-vanishing "effective string tension" (^^ depending on the external 
field (being in deconfinement at zero field). However, this behavior of the Polyakov loop 
does not indicate a restoration of confinement for certain external field fiuxes. At special 
fiux values (for which the internal electric field is quantized according to eq. (^)) both 
Polyakov loop expectation value and cxeff coincide with the values at zero external field. 

The string tension (correctly defined from the correlation function of Polyakov plane- 
plane correlators parallel to the external electric field in spatial direction) is not infiuenced 
by the external electric field and coincides with the zero-field value. 
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The tree level effects on Polyakov loops and Polyakov loop correlation functions are 
absent for external magnetic fields pointing in the time-like direction. 

The dynamics of the Abelian system is different from the behavior of the (3 + 1)D 
non-Abelian gauge theory. The authors of Ref. have found that the external magnetic 
field increases the deconfinement temperature contrary to our results in {2 + l)D compact 
Abelian gauge theory. The reason of this difference may lie in the different behavior of 
the monopoles in the Abelian and non- Abelian gauge theories. The investigation of the 
monopole properties in the non-Abelian gauge theory is underway p4 . 



Acknowledgements 

M. N. Ch. acknowledges a support of Sachsisches Staatsministerium fiir Kunst und Wis- 
senschaft, grant 4-7531.50-04-0361-01/16 and kind hospitality of NTZ and the Institute 
of Theoretical Physics of Leipzig University. Work of M. N. Ch. was partially supported 
by grants RFBR 99-01230a, RFBR 01-02-17456, INTAS 00-00111 and CRDF award RPl- 
2103. E.-M. I. acknowledges the support by the GraduiertenkoUeg Quantenfeldtheorie for 
a working visit to Leipzig. 



References 

A. M. Polyakov, Nucl. Phys. B120 (1977) 429. 
M. Gopfert and G. Mack, Commun. Math. Phys. 82, 545 (1981). 

B. Svetitsky, Phys. Kept. 132 (1986) 1. 
P. D. Coddington et ai, Phys. Lett B175 (1986) 64. 
N. Parga, Phys. Lett. B107 (1981) 442. 

J. M. Kosterlitz and D. J. Thouless, J. Phys. C6 (1973) 1181. 

M. N. Chernodub, E.-M. Ilgenfritz and A. Schiller, |iep-lat/0105021| , to appear in 
Phys. Rev. D. 



M. N. Chernodub, Phys. Rev. D63 (2001) 025003; |hep-th/001lT2l 



B. L. Bakker, M. N. Chernodub and A. 1. Veselov, Phys. Lett. B502 (2001) 338; 
I. I. Kogan, A. Kovner and B. Tekin, JHEP 0103 (2001) 021. 

N. O. Agasian and K. Zarembo, Phys. Rev. D57 (1998) 2475. 

G. Dunne, I. I. Kogan, A. Kovner and B. Tekin, JHEP 0101 (2001) 032. 

T. Schafer, E. V. Shuryak, J. J. M. Verbaarschot, Phys. Rev. D51 (1995) 1267; 

R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Annals Phys. 280 (2000) 35; 

E.-M. Ilgenfritz and E. V. Shuryak, Phys. Lett. B325 (1994) 263. 



21 



[12] N. O. Agasian, Phys. Lett. B488 (2000) 39; 

1. A. Shushpanov and A. V. Smilga, Phys. Lett. B402 (1997) 351; 

A. I. Vainshtein et ai, Sov. J. Nucl. Phys. 39 (1984) 77; 

R. Brower, P. Rossi and C. Tan, Nucl. Phys. B190 (1981) 699. 

[13] R Cea and L. Cosmai, Phys. Rev. D48 (1993) 3364; |hep-lat/01010T7 . 

[14] M. N. Chernodub, F. V. Gubarev and E.-M. Ilgenfritz, Phys. Lett. B424 (1998) 106; 
M. N. Chernodub, F. V. Gubarev, E.-M. Ilgenfritz and A. Schiller, Phys. Lett. B443 
(1998) 244; ihid. B434 (1998) 83. 

[15] K. Kajantie et. al, Nucl. Phys. B544 (1999) 357; 

D. Comelh, D. Grasso, M. Pietroni and A. Riotto, Phys. Lett. B458 (1999) 304; 
V. Skalozub and M. Bordag, Int. J. Mod. Phys. A15 (2000) 349. 

[16] P. Cea and L. Cosmai, Phys. Rev. D48 (1993) 3364. 

[17] P. Dimopoulos, K. Farakos and G. Koutsoumbas, Phys. Rev. D63 (2001) 054507. 

[18] J. Alexandre et ai, |hep-lat/01010Tl| ; 

K. Farakos et ai, Phys. Rev. D61 (2000) 045005; 

V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, Phys. Rev. D52 (1995) 4747; 
Phys. Rev. Lett. 73 (1994) 3499 [Erratum-ibid. 76 (1996) 1005]. 

[19] P. H. Damgaard and U. M. Heller, Nucl. Phys. B309 (1988) 625. 

[20] P. Cea and L. Cosmai, Phys. Rev. D60 (1999) 094506. 

[21] P. Cea and L. Cosmai, Phys. Rev. D43 (1991) 620. 

[22] J. Glimm and A. Jaffe, Comm. Math. Phys. 56 (1977) 195. 

[23] J. Engels, V. K. Mitrjushkin and T. Neuhaus, Nucl. Phys. B440 (1995) 555. 

[24] M. N. Chernodub, E.-M. Ilgenfritz and A. Schiller, in preparation. 

[25] T. A. DeCrand and D. Toussaint, Phys. Rev. D22 (1980) 2478. 

[26] Z. Schram and M. Teper, Phys. Rev. D48 (1993) 2881. 



22 



